NPS-62-80-017PR 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


AN  ADAPTATION  OF  THE  MINIMUF  HF  PROPAGATION 
PREDICTION  PROGRAM  TO  THE  TI-59  CALCULATOR 


Jay  H.  Benson 
Stephen  Jauregui,  Jr. 


Final  Report  for  period  April  1980  -  July  1980 


Approved  for  public  release?  distribution  unlimited. 

Prepared  for:  Naval  Electronic  Systems  Command 
(PME-107-6) 

Washington,  DC  20360 


20091105019 


7  K 

s/ox 

846 


1 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey)  California 


Rear  Admiral  J.J.  Eke  lurid 
Superintendent 


Jack  R.  BorstinS 
Provost 


The  work  reported  herein  was  supported  by  the  Naval 
Electronic  Systems  Command* 

Reproduction  of  all  or  part  of  this  report  is 
authorized ♦ 

This  report  was  prepared  by? 


Lieutenant)  United  States  Navy 


STEPHEN  JAUREGUIr  Jr. 

Associate  Professor  of  Engineering 


Reviewed  by 


Released  byJ 


DONALD  E.  KIRK 
Chairman)  Department  of 
Electrical  Engineering 


WILLIAM  M.  TOLLES 
Dean  of  Research 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Dmtm  Entmrmd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  *  2.  GOVT  ACCESSION  NO. 

NPS -62-80-01 7PR 

3.  RECIPIENT'S  CATALOG  NUMBER 

«.  TITLE  (and  Submit) 

AN  ADAPTATION  OF  THE  MINIMUF  HF  PROPAGATION 
PREDICTION  PROGRAM  TO  THE  TI-59  CALCULATOR 

S.  TYPE  OF  REPORT  8  PERIOD  COVERED 

Final  Report  April-July  1980 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORf*; 

Jay  H,  Benson 

Stephen  Jauregui,  Jr. 

§.  CONTRACT  OR  GRANT  NUMBER^) 

9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 

Monterey,  California  939^-0 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

II  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Electronic  Systems  Command  (PME-107-6) 
Washington,  DC  20360 

IZ.  REPORT  DATE 

1  July  1980 

13.  NUMBER  OF  PAGES 

27 

14  MONITORING  AGENCY  NAME  A  AODrESS(//  dlf format  from  Controlling  0//ic#> 

IS.  SECURITY  CLASS,  (ot  OH.  rmporl) 

UNCLASSIFIED 

is*.  DECLASSI  F|  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thlm  Rmport) 


Approved  for  public  release;  distribution  unlimited# 


17.  DISTRIBUTION  STATEMENT  (of  thm  mbmtract  mntmrmd  In  Block  20,  If  dllfmrmnt  from  Rmporl) 


IB,  supplementary  notes 


19  KEY  WORDS  (Contlnuo  on  rmvmrmm  mldm  If  nmcmmmmry  and  Identify  by  block  numbmr) 

MINIMUF 

HF  PROPAGATION  PREDICTION 
MAXIMUM  USABLE  FREQUENCY 


20.  ABSTRACT  fConNnu#  on  rmvmrmm  mi  da  If  nmcmmmmry  and  Idmntlty  by  block  numbmr) 

An  algorithm  for  predicting  the  Maximum  Usable  Frequency  between  two 
points  on  the  surface  of  the  earth  for  High  Frequency  communications 
is  implemented  on  a  TI-59  hand-held  programmable  calculator. 


DD  ,  j°nM7!  1473  edition  OF  I  NOV  88  IS  OBSOLETE  UNCLASSIFIED 

S/N  0  10  2- 014-  660 1  | 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Dmtm  Mntmrmd) 


TABLE  OF  CONTENTS 


ABSTRACT 


1.  INTRODUCTION  -  - - _______________  3 

2.  OBJECTIVE  - - ___________  3 

3.  PROGRAM  DEVELOPMENT - - - ______  4 

4.  PERFORMANCE  4 


5.  CONCLUSIONS  AND  RECOMMENDATIONS 


APPENDIX  A  BASIC  MINIMUF-3.5  PROGRAM  LISTING  -  -  6 

APPENDIX  B  TI-59  PROGRAM  LISTING  -------  -  9 

APPENDIX  C  PROGRAM  LOGIC  MODIFICATIONS  -----  ig 

APPENDIX  0  PROGRAM  OPERATING  INSTRUCTIONS  -  -  -  20 

APPENDIX  E  TEST  CASE  --------------  22 

APPENDIX  F  TI-59  CUSTOM  MODULES  --------24 

LIST  OF  REFERENCES  - - - - - ------25 

INITIAL  DISTRIBUTION  LIST  -  _____  -  -  -  -  26 


1 


ABSTRACT 


An  algorithm  for  predicting  the  Maximum  Usable  Frequency 
between  two  points  on  the  surface  of  the  earth  for  Hislh 
Freauenca  communications  is  implemented  on  a  TI-59  hand-held 
programmable  calculator. 


1*  INTRODUCTION 


An  algorithm  for  predicting  the  Maximum  Usable 
Frequency  (MUF)  in  the  High  Frequency  radio  band  has  been 
developed  at  the  Naval  Ocean  Systems  Center t  San  Diego* 
California*  by  R.B.  Rose  and  J.N.  Martin  CH*  and  has  been 
implemented  in  both  FORTRAN  IV  and  BASIC  computer  languages 
under  the  name  MINIMUF-3.5.  The  algorithm  has  been 
incorporated  into  the  CLASSIC  F'ROF'HET  propagation  analysis 
system  *  and  has  been  demonstrated  to  provide  a 
field-deployable  capability  for  computing  HF  propagation 
paths  with  micro-computer  resources. 

2.  OBJECTIVE 

The  essential  portions  of  MINIMUF-3.5  are  contained  in 
approximately  80  BASIC  language  statements  *  reproduced  in 
appendix  A*  which  suggests  a  potential  for  being  encoded  in 
hand-held  programmable  calculators.  H.F.  Hite  of  Hughes 
Aircraft  Company  adapted  a  restricted  version  of  MINIMUF-3.0 
for  the  HP-67  calculator  that  computes  MUFs  for  1-hop 
F-laaer  propagation  and  is  thus  limited  to  station 
separations  of  4000  kilometers  or  less.  A  complete 
implementation  of  MINIMUF*  capable  of  2-hop  path  prediction* 
was  desired. 
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3.  PROGRAM  DEVELOPMENT 


A  TI-5?  calculator  program  was  adapted  directly  from 
the  Rose  and  Martin  BASIC  program  listing  with  some  logic 
changes  reouired  to  accomodate  the  memory  constraints  of  the 
calculator*  The  program  listing*  summary  of  logic  changes* 
operating  instructions*  and  test  case  are  contained  in 
appendices  B*  C*  D*  and  E*  respectively*  The  program  listing 
is  annotated  with  cross-references  to  statements  in  the 
BASIC  listing  for  assistance  in  tracing  logic  flow.  The  hulk 
of  the  development  task  was  to  minimize  program  and  data 
storage  reaui rements ♦  Several  iterations  finally  produced 
the  current  version  which  reaui res  800  program  steps  and  20 
data  registers*  the  full  capacity  of  the  normal ly-accessible 
calculator  attributes* 

4.  PERFORMANCE 

The  TI-59  program  was  optimized  for  storage*  that  being 
the  overriding  constraint.  Execution  time  for  single-hop 
predictions  is  approximately  50  seconds.  Were  the  program 
optimized  for  speed*  given  a  larger  storage  capacity* 
execution  time  could  conceivably  approach  40  seconds*  The 
test  case  contained  in  appendix  E  is  the  identical  test  case 
promulgated  with  the  Rose  and  Martin  report*  and  yields  the 
same  results  on  the  TI-59*  Notice*  however*  that  the  MUF 
limit  of  32  MHz  has  been  raised  to  50  MHz*  thereby  providing 
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for  HF  path  prediction  during  periods  of  high  solar  flux 
when  customary  band  limitations  are  exceeded.  The  test  case 
should  be  executed  after  initially  programming  the  TI-59  by 
keystroke  to  ensure  correct  program  entry.  A  copy  of  the 
program  is  available  on  magnetic  cards  by  sending  two  blank 
cards  to  the  authors*  however*  the  user  is  cautioned  that 
magnetic  cards  are  not  guaranteed  to  be  transportabl e  among 
all  TI-59  calculators. 

5.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  capability  for  predicting  Maximum  Usable  Frequencies 
in  the  HF  radio  band  can  be  realized  with  hand-held 
’ calculators  *  and  therefore  can  be  widely  deployed  in  the 
field.  An  even  larger  capability  can  be  realized  by 
utilizing  custom  made  modules  for  the  TI-59.  Each  module 
provides  for  5000  program  steps  in  addition  to  the  normal 
960  maximum  step  capacity  of  the  calculator*  and  could 
provide  for  additional  portions  of  the  CLASSIC  PROPHET 
System  to  be  included  in  calculator  implementations  *  such  as 
the  D-region  absorption  model.  Appendix  F  lists  several 
persons  to  contact  in  regard  to  having  custom  modules 
manufactured  for  the  TI-59. 
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APPENDIX  A 


BASIC  MINIMUF-3 . 5  PROGRAM  LISTING 


The 

BASIC 

language  listing 

of  MINIMUF-3. 5  is 

reprinted  here 

with 

the 

permission  of  R.B 

.  Rose*  and  serves 

as  a  guide  for 

the 

T I -59 

keystroke  listing 

♦ 
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MINIM  til- -0.5  PUOCiK  VM 


1000  REH 

1010  K7-SIN<L1>*SIN<L2>+C0S<L1>*C0S<L2>*C0S<M2-M1> 

1020  G1-ACS<X7  MAX  -1+1.0E-5  MIN  1-1.0E-8) 

1038  X6*1.59*G1 
1040  K6-K6  MAX  1 
1030  K5*l/K6 
1060  J9*100 

1070  FOR  K1-1/<2*K6>  TO  l-l/<2 *X6>  STEP  0.9999-1/X6 
1080  IF  K3-1  THEN  1100 
1090  K3-0.3 
1100  P-SIN<L2> 

1110  Q*C0S<L2> 

1129  A-<SIN<Lt>-P*COS<Gl>>/<Q*SIN<Gl>> 

1139  B*G1*K1 

1140  C-PSCOS(B>+Q*SIN<B>*A 

1150  D*<C0S<B>-C*P>/<Q*SQR<1-C*C>> 

1160  D«ACS<D  MAX  -1+1.0E-5  MIN  1-1.0E-8) 

1170  M0«M2+SGN<SIN<M1-W2>>*D 
1180  IF  W8«>0  THEN  1200 
1196  M8-H0+P1 
1208  IF  U0<P1  THEN  1220 
1210  U0-H0-P1 

1220  L0*P0-ACS<C  MAX  -1+1.0E-5  MIN  1-1.0E-8) 

1238  Y1»0.0172*<10+<M0-1>*30.4+D6> 

1248  Y2»0.409*COS<Y1> 

1250  X8*3 . 823 M0+ 1 2+0 . 13*<SIN<Yl>+i.2*SIN<2*Yl>> 

1268  X8*X8-123< l+SGN<K8-24> ><SGN<ABS<K8-24>  > 

1270  IF  COS <L0+Y2> >-8.26  THEN  1350 

1288  X9*0 

1290  G0-0 

1380  M9-2.53G13X3 

1310  M9-M9  MIH  P0 

1320  M9*SIH(M9> 

1338  M9* 1 *2 . 33M93SQR <  M9  > 

1348  GO  TO  1590 

1330  X9-<-0. 26+SIN<Y2>*SIN<L0>  >/<COS(Y2>*COS(L0>+l . 8E-3> 

1360  X9-12-ATN<X9/SQR<ABS<1-X9*X9>> >*7.639437 

1370  T-K8-K9/2+12*<l-SGN<K8-K9x2>>*SGN<ABS<K8-K9/2>> 

1380  T 4*K8+K9/'2- 1 2$  <  1  ♦SGN < X8+K9/2-24  >  >  &SGN  <  ABS  <  X8+X9/2-24  )  > 
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1390  C0*ABS<COS<L0+Y2>> 

1400  T9*9. 7*C0t9. 6 
1410  IF  T9>0. 1  THEN  1430 
1420  T9=0. 1 
1430  M9-2. 5*G1*K5 
1440  H9*H9  HIM  P0 
1450  H9*SIN<H9> 

1460  H9*  1  +2 . 5*H9*SQR  <  H9  > 

1470  IF  T4<T  THEN  1500 

1480  IF  <T5-T)*<T4-T5)>0  THEN  1510 

1490  GO  TO  1640 

1500  IF  <T5-T4>*<T-T5>>0  THEN  1640 

1510  T6ST5+12$< 1+SGN<T-T5> )$SGN< ABS<T-T5) ) 

1529  G9«PI*<T6-T>/K9 

1530  G8*PItT9/K9 
1540  U*<T-T6>/T9 

1550  G0*C0*<SIN<G9>+G8*<EXP<U)-COS(G9> ) >/( 1+G8*G8> 

1568  G7*C0*<G8*<EXP<-K9/T9>+1 ) >*EXP<  <K9-24>/2>/< 1+G8*G8> 
1570  IF  G0«>G7  THEN  1590 
1580  G0*G7 

1590  G2~<  1  +S9/250 > TH9TSQR < 6+581SQR < G0 )  ) 

1600  G2«G2*< 1-0. 1*EXP< <K9-24>/3>  > 

1610  G2-G2*<1+<1-SGN<L1)*SGN<L2>>*0. 1) 

1629  G2-G2*< 1-0.  1*< 1+SGN< ABS<SIN<L0>  >-COS<L0>  > ) ) 

1630  GO  TO  1700 

1640  T 6*T5+ 1 2$ < 1 +SGH < T 4-T5 ) >*SGN<ABS<T4-T5>  > 

1650  G8*PI*T9/K9 
1668  U-<T4-T6>/2 
1670  U1—K9/T9 

1680  G0*C0t ( G8t < EXP <  U 1  )  +  l  )  >*EXP<U>/<  1+G8*G8> 

1690  GO  TO  1590 

1780  IF  G2>J9  THEN  1720 

1710  J9*G2 

1720  NEXT  K1 

1730  J9-J9  MAX  2  HIN  32 

1740  RETURN 
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APPENDIX  B 


TI-59  PROGRAM  LISTING 


The  keystroke  listing  of  the  TI-59  program  implementation  of 
MINIMUF-3.5  follows.  Segments  of  TI-59  code  are 
cross-referenced  to  the  BASIC  prodram  listing. 
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APPENDIX  C 


PROGRAM  LOGIC  MODIFICATIONS 


The  following  program  logic  changes  were  made  to  the  BASIC 
version  of  MINIMUF-3.5  to  minimize  program  storage 
reoui rements J 

(1)  The  computation  of  M9  at  lines  1300-1330  and  1430-1460 
is  independent  of  the  main  loop  and  was  moved  to  the 
beginning  of  the  program*  immediately  following  the 
calculation  of  G1  and  K5  * 

(2)  The  loop  index  computation  at  line  1070  was  separated 
into  a  computation  of  a  new  variable*  KH0P*  whi-ch  is  either 
.1  or  2*  and  a  follow-on  computation  of  K1  from  KH0P  as  the 
first  item  within  the  loop* 

(3)  The  computation  of  G8  at  lines  1530  and  1650  is  common 
to  both  legs  of  the  preceedinS  test  for  < T5-T4 ) < T-T5 ) *  and 
h3s  been  moved  ahead  of  that  test. 

(4)  The  test  at  line  1470  serves  to  reverse  the  sense  of  the 
following  tests  at  lines  1480  and  1500.  This  logic  has  been 
combined  into  a  single  test  on  < T5-T4 ) < T-T5 )  with  a 
weighting  factor  of  1  or  -1  to  reverse  the  sense  of  that 
test . 
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(5)  The  computation  at  line  1360  was  simplified  to  the 


arcsine  by  way  of  trigonometric  identity. 

(6)  The  MUF  limit  established  at  line  1730  was  changed  to  50 
vice  32  to  accomodate  high  solar  flux  densities. 

(7)  Throughout  the  program?  explicit  calculation  and  storage 
of  variables  that  are  used  only  once  in  following  statements 
has  been  eliminated  to  conserve  on  data  register 
reou i rements . 
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APPENDIX  D 


PROGRAM  OPERATING  INSTRUCTIONS 


The  following  instructions  must  he  followed  to  operate  the 
TI-59  version  of  MINIMUF-3.5J 

(1)  Repartition  the  calculator  for  800  program  steps  and  20 
data  registers  by  pressing  2/0P/17. 

(2)  Load  the  4  memory  banks  from  program  cards  (2)  or  by 
entering  program  keystrokes  with  the  calculator  in  the  LRN 
mode  ♦ 

(3)  Enter  input  data  as  follows? 

(a)  Transmitter  North  Latitude  in  decimal  degrees  in 
R12  (range  -90  degrees  to  90  degrees) 

(b)  Transmitter  West  LonStitude  in  decimal  degrees  in 
R13  (range  0  to  360  degrees) 

(c)  Receiver  North  latitude  in  decimal  degrees  in  R14 
(range  -90  to  90  degrees) 

(d)  Receiver  West  Longtitude  in  decimal  degrees  in  R15 
(range  0  to  360  degrees) 
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(e)  Month  in  R16  ( rande  1  to  12) 

(f)  Day  in  R17  ( rande  1  to  31) 

(S)  Sunspot  Number  in  R19  ( ranSe  -  positive  number) 

(4)  Enter  Time  in  X-reSister  ( ranSe  0  to  24  decimal  hours) 

(5)  Press  A 

(6)  If  a  printer  is  attached*  the  input  time  is  echoed* 
followed  by  the  answer  (MUF)  after  50  seconds  for  sinSle-hop 
and  100  seconds  for  two-hop  predictions* 

(7)  The  answer  is  displayed  in  the  X-reSister  when  the 
program  halts* 

(8)  All  input  Quantities  remain  undisturbed  in  R12  thru  R19 
(time  is  stored  by  the  program  in  R18>*  Another  MUF  for  a 
different  time  can  be  computed  directly  by  repeating  steps 
(4)  and  (5)* 
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APPENDIX  E 


TEST  CASE 


The  test  case  that  follows  was  provided  in  the  NOSC  report 
on  MINIMUF-3 . 5 .  The  actual  printer  listing  of  input  and 
program  output  is  included. 
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APPENDIX  F 


TI-59  CUSTOM  MODULES 


The  following  points  of  contact  are  provided  for  the  design 
and  manufacture  of  TI-59  custom  modules*  Both  Texas 
Instruments  and  Horizons  Technology  (contracting  with  TI) 
provide  software*  emulator*  and  consulting  support  for  the 
design  and  production  of  custom  modules*  Costs  Quoted  at 
this  time  are  approximately  $12*000  for  a  minimum  order  of 
250  modules*  This  includes  emulator  support  and  manufacture* 
but  does  not  include  software  consulting*  Cost  reduction 
through  Quantity  is  available* 


Fred  Wilke 

TI-59  Custom  Module  Division 
Texas  Instruments 
Lubbock*  TX 
(806)  741-3240 

Robert  Kruger 

Horizons  Technology*  Inc* 

7830  Clairemont  Mes3  Boulevard 
San  Diego*  CA  92111 
(714)  292-8331 
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